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Dedicated to Professor Jon Borwein on the occasion of his sixtieth birthday. 

Abstract. We give an example of a formula involving the sine function that holds for every 
N = 0,1,2, .. ., up to about ioi''2832732i65^ then fails for all larger N. We give another example 
^— I that begins to fail after about iV ~ exp(exp(exp(exp(exp(exp(e)))))). This number is larger 

than the Skewes numbers. 
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^ 1. Introduction 

T— I 

^ In a 1992 paper [6], Jon and Peter Borwein give examples of formulas that are "almost" true: 

• ^ that is, they are correct to anywhere from thousands to over 42 billion decimal places, but are 

^ not actually true. Here, we'll do something similar. We'll give examples of a formula involving 

d the sine function that holds for a ridiculously large number of values of A = 0, 1, 2, . . . before 

it begins to fail. 

The sine function is defined as: sinc(O) = 1, and sinc(a;) = sin(a;)/a; if x 7^ 0. 
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The sine funetion has many interesting properties. For example, there are the curious iden- 
tities 



oo 



sinc(n) = sinc(n 

and 



n=l n=l 



sinc(a;) dx = sinc(x)^ dx 
Jo 
Both sums equal 7r/2 — 1/2; both integrals equal 7r/2; see [2]. 

David and Jon Borwein and Bernard Mares [5] showed how to evaluate certain integrals 
involving products of sine functions. A more recent paper [2] shows a connection between sums 
and integrals of products of sine functions. That paper explains the curious identities above, 
and further explains why the integrals are 1/2 more than the sums. 

Let > and Oq, ai, a2, . . . be positive numbers. Theorem 1 of [2] states, in part, that 

N N 



Y\_smc{akx) dx (1.1) 



smc[akn) = 

n=l k=0 " ^ k=0 

provided that 

TV 

^afc<27r. (1.2) 

k=0 

If = 0, the condition required for equality in (1.1) is simply that Qq < 2tt. 

If is a slowly divergent series, then the sum in (1.2) will exceed 27r only when A^ exceeds 
some very large number Aq- Therefore, equation (1.1) will be true for the many cases A^ = 
0, 1, 2, . . . A^O; and will fail for all A^ > Aq. We will show how to construct interesting examples 
with arbitrarily large Nq. We do this by finding infinite series that diverge very, very slowly. 

For example, when the at are the reciprocals of primes of the form lOn + 9, then we can 
estimate that A^q ~ I0i02832732i65_ 

We could make estimates in a similar way if the sum in (1.2) is a convergent series whose sum 
exceeds 27r, but we will not consider that case here. Our results follow from standard theorems 
in number theory and numerical analysis. 

We should also make clear that, even if we are able to get only a rough estimate for the value 
of A'o such that equation (1.1) fails for N > Nq, this cutoff is, nevertheless, well-defined. 

2. An Example With Odd Denominators 

Example 1 (a) in [2] illustrates Theorem 1 of that paper using the sequence = l/{2k + 1) 
for A; > 0. In this case, one can calculate that 

N 
k=0 

does not exceed 27r until A^ > 40249. Therefore, equation (1.1) holds for every A^ = 0, 1, 2, . . ., 
40248, but fails for A^ > 40249. Interestingly, for this example, Crandall has shown [8, p. 24] 
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that, at = 40249, the left side minus the right side of (1) is positive but less than 10 226576_ 
So, even if one did use a computer to check whether 



^ oo N 



(2F?T) = rn;'-fe)''" 



sine 

n=l fc=0 " ■ -/ -'u i^^Q 

all the way up to = 40249, the left and right sides would appear to be the same at = 40249 
unless one had the foresight to perform the calculation to over 225,000 decimal places! 

We all know that it is dangerous to rely on a formula merely because it is true for a few test 
cases, or a few hundred, or even a few thousand. However, everyone has done this at one time 
or another. This example shows how truly dangerous that practice can be! 

3. More Examples With Denominators in Arithmetic Progressions 

There's nothing special about the "2" in the denominators 2A; + 1 in the example above. 
What happens with a^. = l/{mk + 1) for other values of m? 

As before, (1.1) will hold as long as (1.2) holds. For m = 1, 2, 3, and 4, one can add up terms 
of the series until the sum exceeds 27i. However, as m increases, the number of terms required 
to make the partial sum exceed 27r increases quite rapidly. For larger values of m, we can use 
the Euler-MacLaurin summation formula to accurately estimate the sum. This enables us to 
find the exact number of terms required to make the partial sum exceed 27r. The following 
version of the Euler-MacLaurin summation formula is based on taking a = 0, 6 = n in [7, p. 
309]: 



E 

A:=0 



fik)= I /(x)rfx+^(/(0) + /(M)) (3.1) 







t (2J)! 



Bk is the k Bernoulli number. The remainder R is bounded by 

where xq is some number between and M, inclusive. Note that this gives us an approximation 
to the sum of the first M + 1 terms (not M terms) of the series. 

A naive approach would be to apply Euler-Maclaurin summation to the function f{x) = 
l/{mx + 1). Then the derivative of order (2J + 2) of f{x) will be of the form C /{mx + 1)^"'+'^. 
The error term R achieves its maximum value over [0, M] at a; = 0. This maximum value is 
not necessarily small. Often, it is too large for the formula to be useful. 

For example, when we apply the Euler-Maclaurin formula with J = 1 to f{x) = l/{2x + 1), 
we get 

*^ 2 1 11 

/(^) - o + WTWTTT-^ - WWTT—^ + o log(2M + 1) 



3 2(2M+1) 6(2M+1)2 2 
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with an error term of 

im < 



15(2x + 1)5' 

However, watch what happens if we separately compute the sum of the first 100 terms of the 



series 

99 



s = V — ^ = 1 + 1 + 1 + . 3.28434218930163434565, 

A:=0 

and then apply the Euler-Maclaurin formula to the "tail", f{x) = l/{2x + 201). We get 

M 



fc=0 



2k + 201 

302 1 1 1, . X log(201) 



121203 2(2M + 201) 6(2M + 201)^ 2 ' 2 ' 

The error term is now 

, , 8M , , 

R < — -. 3.2 

' ' - 15(2a; + 201)^5 ^ ' 

For a given M, this error term is much smaller at x = than was the previous expression for 
R. 

So, in order to estimate the (largest) number of terms for which 



fc>0 



we compute the value of M such that 



.„., + ^log(2M + 201)-™ 



121203 2(2M + 201) 6(2M + 201)^ 2 '2 

equals 27r - s ~ 2.99884311787795213128. The solution of this equation is M ~ 40148.81104. 
But to be sure that M = 40148 is really the value we want, we must verify that: 

• (1) using M = 40148, the estimated sum < 27i 

• (2) the estimated sum plus the next term > 27t 

• (3) the error term is less than 1 / (last term) 

• (4) the error term is less than 2it - (estimated sum) 

Evaluated at M = 40148 and x = 0, error term (3.2) is about 6.52653 ■ 10^^. This is less than 
any term in the series near this value of M. Moreover, the sum of the initial terms plus the 
estimated value of the next M + 1 terms, differs from 2tt by about 1.008 ■ 10~^. 
Therefore, it is safe to conclude that 

99 40148 40248 ^ 

' < 27r 



5^ 2A; + 1 ^ ^ 2A; + 201 2k ^ 

k=0 k=0 k=0 
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and 

99 ^ 40149 ^ 40249 ^ 

5^ 2k + 1 ^ ^ 2k + 201 ~ ^ 2k + 1 ^ 

fc=0 fc=0 fc=0 

Therefore, the largest value of M such that 

M ^ 

< 2tt 



E 



2k + 1 
is M = 40248. 

For 171 = 1 through m = 20, we compute the largest value of M such that 

M 



E 



1 

< 27r. 



mk + 1 

fc=0 

We use the Euler-Maclaurin formula as in the previous example. For larger m, we need higher- 
order approximations in order to compute M with sufficient accuracy and to be assured that 
the error term is small enough. For each m, Table 1 shows the value of M, along with K, the 
number of initial terms, and the value of J that was used in the Euler-Maclaurin formula. With 
a sufficient number of initial terms and a modest value of J, we can keep the error small. 

In all cases, the "tail" function f{x), to which we apply the Euler-Maclaurin formula, is given 
by f{x) = l/{mx + Km + 1). 

(We do not claim that the choice of either 7^ or J is minimal or, in any sense, optimal. We 
could use fewer initial terms, but we would then need to increase J in order to keep the error 
term small enough. Making J larger increases the complexity of the equation involving M.) 

As a check on this procedure, when m = 1, this technique successfully finds the 44-digit value 
of M obtained in [4] for which 

M 

y— 

^ k + l 

exceeds 100. See also [12]. This value can be calculated, for example, using K = 10000 and 
J = 10. We can also verify the first four values of M in the table by direct summation. 

As a further check, if we double the number of initial terms and increase J by 1 or 2, then 
we get the same values of M shown in Table 1. 

We conclude that, for each m = 1, 2, . . . , 10, the equation 

^ oo A'" f 71 \ ^ ( \ 

n=l fc=0 ^ ^ •'^ k=0 ^ ^ 

holds precisely for those iV = 0, 1, 2, . . . M, where M is the corresponding value in Table 1, 
and that it fails for larger iV. Notice that, for m = 20, M exceeds 10^^. So, if we were testing 
"random" values of A^, looking for a value that made (1.1) fail, we would be unlikely to find 
one. 

The exact value of M corresponding to m = 100 can be obtained from the Euler-Maclaurin 
formula by taking K = 50000 and J = 60. The result is the 230-digit value of M shown below. 
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m 


M 


K 


J 


1 


299 


100 


1 


2 


40248 


100 


1 


3 


6699356 


1000 


1 


4 


1199640415 


1000 


2 


5 


222209538173 


1000 


3 


6 


41928392459412 


1000 


3 


7 


8000791810720605 


1000 


4 


8 


1537961933686185453 


1000 


5 


9 


297136851932007766218 


1000 


6 


10 


57616626381701142703593 


1000 


7 


11 


11202463675353183172586907 


1000 


8 


12 


2182608774487516995740392959 


1000 


8 


13 


425930131275278060684396950683 


1000 


10 


14 


83225800344072649528263059652618 


1000 


10 


15 


16279019516889202909861702224716180 


1000 


11 


16 


3186898150182578894413446451622161442 


1000 


12 


17 


624331345650550634164994069452043597341 


1000 


13 


18 


122382651928233262387099042295616064808177 


1000 


14 


19 


24001760343280992647777613927571508451532532 


1000 


15 


20 


4709265577657827035628502288018792360631413283 


1000 


16 



Table 1. The largest M for which ^^^^^ < 2n 



M =15930636153764656093549951961696713197434975028940 (3.3) 
85877192998763567162101035983381719598376913882972 
95285352168437589967676947222915769714257521188927 
15116548003599042566741587106668007049302125094673 
665769807765071841758755530945. 

Therefore, we know that, for = 1/(100A; + 1), equation (1.1) holds for N < M, and fails 
for N> M. 

4. An Example With Primes 

Example (2) in [2] is even more striking. Here we pick to be the reciprocals of the primes, 
with ao = 1/2, ai = 1/3, etc. The sum of the reciprocals of the primes 



fc>0 k>0 ^'^ 



FUN WITH VERY LARGE NUMBERS 



7 



diverges. Then (1.1) becomes 

TV 

2+ V 



^ oo TV ^ \ ^ / \ 



n=l A:=0 ^ "^^ fc=0 

which will hold only as long as 

^ 1 

So, how long is that? From analytic number theory (see, e.g, [9, p. 35] or [1, p. 156] [AP], 
we know that 



log logx + B + o(l) 



where B ~ .26149 ... is Mertens' constant. For computational purposes, we will simply drop 
the o(l) term. In order to make this sum reach 2tt, we must have 

X ~ exp(exp(27r - 5)), (4.2) 

so that X ~ 10^''^. Then, the where (4.1) ceases to hold is the number of primes up to that 
X, which, by the Prime Number Theorem [9, p. 10], is roughly x/ logx ~ 10^^^. 

Without assuming the Riemann hypothesis, Crandall [8, Theorem 11] has shown that, at the 
first for which (4.1) fails, the left side minus the right side is positive but less than 10^*^^° \ 
Now, if both sides of (4.1) were on the order of 10"^^°^^^^ then this small difference would not 
be very interesting. However, neither side of (4.1) is tiny. Crandall calculates [8, Theorem 4] 
that the right side exceeds .6867r for all A^. 



5. Some Examples With Primes in Arithmetic Progression 

In the previous section, we saw that equation (4.1) holds "only" up to A^ = A'q — 10^"^^, and 
then fails for larger A^. In this section, we will give examples where the corresponding A"o is 
much larger. To do this, we will use a subset of the primes, the sum of whose reciprocals still 
diverges, but much more slowly than the sum of the reciprocals of all the primes. Specifically, 
we will use the set of primes in an arithmetic progression. 

For positive integers r and s, let (r, s) be the greatest common divisor of r and s. Also, 
(/)(m) denotes Euler's phi function, which is the number of positive integers k < m such that 
{k,m) = 1. 

Dirichlet's theorem states that, if (a, g) = 1, then there are infinitely many primes in the 
arithmetic progression qn + a, where n = 0, 1, 2, ... . The number of primes < a; in such an 
arithmetic progression is asymptotic to 



0(g) log logo; 

We also know that [1, p. 156], if x > 2, 
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E 



1 '°S'°S%A/fe..)+OI^|, (5.2) 



p=a mod g 

where M{q, a) is a constant that depends on q and a. The term inside the 0{...) will be small 
if X is large. However, the constant factor implied by the 0{...) is of unknown size, so we will 
not establish rigorous bounds in what follows. Our goal is only to obtain rough estimates, and 
we will assume that the 0{...) term is small enough to be neglected. 

Given q and a with (a, q) = 1, if the sum up to x equals approximately 27r, then equation 
(5.2) tells us that 

log log X ~ 0(g) • {2tt — M{q, a)), 

so that 

X ~ exp(exp( 0(g) • {2tt — M{q, a)) )). (5.3) 

Compare (5.3) with (4.2). Notice that if g > 3, then 0(g) > 2, which will have the effect of 
multiplying the "top level" exponent by some integer which is at least 2. This means that (5.3) 
will produce much larger values of x (and A^^o) than we got from (4.2). Let's look at several 
examples where M(g, a) is accurately known, so that we can make reasonable estimates. 

Example A. Consider the primes in the arithmetic progression 3n + 1, that is, those primes 
that are = 1 (mod 3). Our sequence will be the reciprocals of these primes: Oq = 1/7, 
Oi = 1/13, a2 = 1/19 fls = 1/31, and so on. The sum of the diverges very slowly. Here, g = 3 
and a = 1. Languasco and Zaccagnini [11] computed the values of many M{q,a) to over 100 
decimals. We will use their value of M(3, 1) rounded to 10 decimals: M(3, 1) ^ —0.3568904795. 
(Using more decimals gives the same final result). Then, since 0(g) = 2, 

log log X = 0(g) ■ (27r - M(g,a)) ~ 2 ■ (27r + 0.3568904795) ~ 13.2801515734 

so 

X ~ exp(exp(13.2801515734)) ~ exp(585459.08163) ~ 4.45176353778 ■ lO^^^aei, 

We now use Equation (5.1) to estimate A^, the number of primes = 1 (mod 3) that are < x. 
We get 

X 4.45176353778-10™! ,^054255 

N ~ — — ~ ~ 3.8 ■ 10254255_ 

0(g) log X 2-585459.08163 

So, for this example, where the in equation (1.1) are the reciprocals of the primes = 
1 (mod 3), equation (1.1) holds for < A^q — 10^^^^^^ and fails for larger A^. 

This number is so much larger than the value we obtained in Section 4, primarily because 
of the factor of 0(g) = 2 in equation (5.3). We were also helped by the fact that M(3, 1) is 
negative. We should also emphasize that we make no rigorous claim about this value of A^ 
because we ignored the 0(...) term in equation (5.2). 

]^g254255 jg large, but we can do better! We'll simply choose an arithmetic progression with a 
value of g that has a larger 0(g). For our next examples, we'll look at the four residue classes 
(mod 10), that is, the primes in the four arithmetic progressions lOn + 1, lOn + 3, lOn + 7, and 
lOn + 9. 
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Example B. Consider the primes in the arithmetic progression lOn + 1. Our sequence will 
be the reciprocals of these primes: Oq = 1/11, oi = 1/31, 02 = 1/41, and so on. Here, g = 10 
and a = 1. Also, (f){q) = 4; observe what a drastic effect this will have on the calculations 
below, compared to Example A. Again, we will use the value of M(10, 1) from [11], rounded to 
16 decimals: M(10, 1) ~ -0.2088344774302376. 

log log X = 0(g) ■ (27r - M(g, a)) ~ 4 • (27r + 0.2088344774302376) ~ 25.9680791384 

so 

X ~ exp(exp(25.9680791384)) ~ exp(l. 89580417544 ■ 10^^) ~ 2.64164832039 • io823337292i6_ 

We now use Equation (5.1) estimate A^, the number of primes = 1 (mod 10) that are < x. 
We get 

X 2.64164832039 ■ 1082333729216 S2ss^72q2n4 

]\[ ^ ^ ~ 3 48 ■ xo 

0(g) log X 4 ■ 1.89580417544 ■ 1011 

So, for this example, where the in equation (1.1) are the reciprocals of the primes in the 
arithmetic progression lOn + 1, equation (1.1) holds until somewhere around A^ ~ ^O*^^'^^''^^^'^^ 
and fails for larger A^. 

Example C. Now take g = 10, a = 3. M(10,3) ~ 0.1386504057476469 and 0(g) = 4. Then 
loglogx = 0(g) ■ (27r - M(g, a)) ~ 4 ■ (27r + 0.1386504057476469) ~ 24.5781396057 

so 

X ~ exp(exp(24.5781396057)) ~ exp(4.72226555917 ■ 10^°) ~ 1.89595583512 ■ io20508538744_ 
Then 

^ ^ ^ 1-89595583512 ■ 10 ^ ^ -|^g20508538733 

0(g) log X 4- 4.72226555917- 1010 
With the lOn + 1 and lOn + 3 sequences, we get estimates near 10^2333729204 ^^^^^ ]^g20508538733^ 
respectively. The calculations were similar. The difference in results arises from the differing 
value of M(g, a). 

Example D. With the arithmetic progression lOn + 7, the first few values are: ao = 1/7, 
ai = 1/17, a2 = 1/37, and ag = 1/47. Rounded to 16 decimals, M( 10, 7) = -0.1039035249178728. 
Carrying out calculations similar to those above, we get 

X ~ 1.38984773649 - io54ii2058088 



and 



2 8 - 10^"' 11 2058076 



iVQ 

Again, the fact that M(10, 7) < helped make the final result, A'q, somewhat larger than that 
for the lOn + 3 progression. 

Example E. With the arithmetic progression lOn + 9, the first few values are: ao = 1/19, 
ai = 1/29, as = 1/59, and ag = 1/49. Rounded to 16 decimals, M(10, 9) = -0.2644151905518937. 
Carrying out calculations similar to those above, we get 

loglogx = 0(g) - (27r - M(g, a)) ~ 4 - (27r + 0.2644151905518937) ~ 26.1904019909 
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SO 

X ~ exp(exp(26. 1904019909)) ~ exp(2.36781116183 ■ 10^^) ~ 9.98876322671 ■ io™2732i76_ 
Then 



X 9.98876322671 ■ 10102832732176 

'0 ~' 



iVn ~ - ~ ^ ^" ~ 1.05 ■ 10l02832732165_ 



0(g) logx 4- 2.36781116183- 10" 
So, for this example, where the ak in equation (1.1) are the reciprocals of the primes in the 
arithmetic progression lOn + 9, equation (1.1) holds for N < Nq c:^ ]^q102832732165 fg^jjg 
larger A^. 
Example F. 

Let's try one final example that has an even larger value of 0(g): q = 100, for which 0(g) = 40. 
We'll take a = 1. The first five primes in the sequence lOOn + 1 are 101, 401, 601, 701, and 
1201. The first few values of are Oq = 1/101, ai = 1/401, 02 = 1/601, and 03 = 1/701. 
Languasco and Zaccagnini [11] computed M(100, 1) to 104 decimals, but, as we shall see, we 
end up with numbers on the order of 10^°^, so 104 decimals appears to be not quite enough. 
Professor Languasco has kindly calculated and provided the following value, accurate to 136 
decimals: 

M(100, 1) ~ -0.0327328506433100964865591320930048072116438944230 
5808121239698784116683056664327790581593738706166 
32469149389219354796589435060666487892. 

Using this more accurate value, and computing x and No as above, we get: 

loglogx = 0(g) ■ (27r - M(g, a)) ~ 40 ■ (27r + 0.0327328506...) ~ 252.6367263129... 

so 

X ~ exp(exp(252.6367263129...)) ~ exp(5.2328244314... ■ 10^°^) ~ 9.1592327310 ■ io22725...82928^ 
where the last exponent on the right has 110 digits. Then 

X 9.1592327310 ■ lO^^''^^'"^^^^'^ ^ ^ _ j^g22725.. .82817 



0(g) logx 40 ■ 5.2328244314... ■ lO^o^ 

where, again, the last exponent on the right has 110 digits. In case the reader is curious, the 
exact value of this integer is 

2272586775359001684288392849910387559794317395514706629 
6853514124083426515979578332298510630142796585419982817. 

Since this is merely an approximation that came from ignoring the 0{...) term, the reader is 
advised against taking this too seriously. An integer with 110 digits is at least lO^''^, so we can 
write 

No ^ 10^°"^ 

Compare this with the value of M given in equation (3.3). There, M was the number of 
terms whose denominators are the arithmetic progression lOOn + 1, and for which the sum 
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remained less than 27r. That M had 230 digits. When the terms are restricted to primes in 
the arithmetic progression 100n + 1, the corresponding number of terms has about 10^''^ digits. 
This is not entirely unexpected: the partial sums of the harmonic series increase as a log; the 
partial sums over primes increase as a log log. 

6. Estimating the Sum of 1/p Over Primes in an Arithmetic Progression 

Before generating more examples with even larger values of Nq, we will discuss equation (5.2) 
in more detail. Equation (5.2) gives an estimate of, but not rigorous bounds, for the size of 

p<x 
p=a mod q 

The reasons for this are twofold. First, the error term is 0(1/ log x), which means that the 
error is bounded by some multiple of 1/ logx. However, we do not know what that multiple is. 
Second, for large q, it is hard to compute M{q, a), so we do not know how large or small M{q, a) 
can be. For q with 3 < g < 300, Languasco and Zaccagnini [11] use Dirichlet L-functions to 
compute M[q, a) 20 decimals for all a with 1 < a < g and (g, a) = 1. All of these numbers are 
available from Languasco's web page; see [11]. 

In this range of q, the largest M value is M(269, 2) ~ .49776. 

For larger q, one can use the approximation 

N loglogX ^1 

p=a mod q 

to get a non-rigorous estimate of M{q,a). For 3 < g < 300, and summing up to x = 10^, 
approximation (6.1) gives values of M{q,a) that agree to 4 or 5 decimals with those in file 
"matricesM.txt" on Languasco's web page. 

For odd g < 10000, M(g, 1) appears to be slightly negative; for example, M(9999, 1) ~ 
-.0004. Also, M(g, 2) appears to be slightly less than .5; for example, M(9999, 2) ~ .49959. In 
fact, this limited data for g < 10000 suggests that M(g, 1) approaches 0, and M(g, 2) approaches 
1/2, as g approaches oo. 

In fact, this is what happens. The paper [11] by Languasco and Zaccagnini uses a result by 
Karl K. Norton to show that, if 1 < a < g and (g, a) = 1, then as g approaches oo, M(g, a) 
approaches 1 /a if a is prime, and approaches otherwise. 

7. Surpassing the Skewes Numbers 
The logarithmic integral of x, li{x), is defined as 

\i(x) = lim / — -— dt + - — -— dt 

^ > e^O \J, log(t) A+. log(t) 

li(x) is a good approximation to tt{x), the number of primes < x. For every value of x for 
which both tt{x) and li(x) have been computed, we observe that 

it{x) < h(x). 
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Nevertheless, Littlewood proved in 1914 that the difference 7i{x) — li{x) changes sign infinitely 
often. In 1933, Skewes [13] proved, assuming the Riemann Hypothesis, that there is an x less 
than the very large number 

79 1 n34 

such that 7r(x) > h(x). 

Then, in 1955, Skewes [14] proved, this time without assuming the Riemann Hypothesis, that 
there is an x less than the much larger number 



such that 7i{x) > li(x). 

More recently. Bays and Hudson [3] use accurate values of the first one million pairs of 
complex zeros of the Riemann zeta function to show that there are values of x near 1.39822 x 
10^^^ such that ^{x) > li{x). In any case, in spite of numerical evidence seemingly to the 
contrary, we do know that there are values of x such that 7r(a;) > li(x). 

These Skewes numbers are much larger than most numbers that are in common use in 
mathematics. 

Here we will present an example where equation (4.1) holds for = 0, 1, 2, . . . A^^o, and fails 
for N > No, where A'^o is much greater than 5*2. The key idea is to use the methods of the 
previous section, but with q such that 0(g) is very large. 

Now let P be the largest known prime, which as of this writing is 2'^^^^^^°^ — 1. This prime 
has 12,978,189 digits, so P > 10^° . We will take our arithmetic progression be the set of 
primes = 1 (mod P). (Since P is odd and we are interested only in primes in this arithmetic 
progression, this is equivalent to the set of primes in the progression 2P + 1). Since P is 
prime, we have 4>{P) = P — 1. Equation (5.3) applies here. The calculations and Norton's 
result mentioned in the previous section suggest that M(g, 1) may be close to 0. We will use 
M(g, 1) = in our calculations. 

With P being the largest known prime, we know from Dirichlet's theorem that the arithmetic 
progression Pn + 1 contains infinitely many primes. We don't know a single example of such a 
prime, but we do know that there are an infinite number of them! 

Let the reciprocals of those primes be our a^. The value of x for which the sum 



surpasses 27r is approximately 

X ~ exp(exp(0(P) • (27r - M(P, 1)) )) ~ exp(exp(27r(P - 1))) > exp(exp(7rP)). 

The last inequality above is not merely an approximation; it holds for P > 2. Therefore, the 
number of primes in the arithmetic progression Pn + 1 up to this x is about 





p=l mod P 



X 



exp(exp(7rP)) 



exp(exp(7rP) — vrP) 



0(P) logx 



(P - 1) exp(7rP) 



P- 1 
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Regarding the numerator, it is easy to convince oneself that if P is large, then e^^ — ttP > e^^ 
(in fact, this holds if p is more than about .391). Therefore, 

exp(exp(7rP) — ttP) exp(exp(2P)) 
° ~ P-1 ^ P-1 

If y > log 2, then > 2 and. Multiplying each side by e^, we get e^^ > 2e^. Taking e to each 
side again, we get 

exp(exp(2y)) > exp(2exp(y)) = exp(exp(y) + exp(y)) = exp(exp(y)) • exp(exp(y)). 
So, if P > 1, 

exp(cxp(2P)) exp(cxp(P)) r r-nw i /m\ 

p_-y = P-\ ' exp(exp(P)) > exp(exp(P)). 
Therefore, we have the very rough estimate, 

A^o ^ exp(exp(P)). (7.1) 

This all works because the factor of P — 1 in the denominator is tiny compared to exp(P) 
and is even smaller when compared to the numerator. Likewise, exp(P) is small compared to 
exp(exp(P)). Our very rough approximation becomes 

A^o ~ exp(exp(P)) ~ exp(exp(10^° )) = ■ (7.2) 

A simple calculation shows that logloglog(P) > 2.845 > e, so that P > e^" . So, we can 
write A^o as 

A'o — exp(exp(P)) ^ 

How large is this A'o? Aq and 5*2, the larger of the two Skewes numbers, are far too large to 
calculate with directly, so we must use logarithms. In fact, we must use logarithms of logarithms 
to bring the numbers within range of most computers. Comparing loglog A^o and log log 5*2, we 
get 

loglogA^o ~ 10^°', 

and 

loglog^s ^ e^"™' ~ 7.6 ■ 10^63 < iQiooo ^ ^qIob^ 

so A'^o is larger than 5*2. How much larger? Because 7.705 ~ e^ °^^^, we can write 5*2 as a tower 
of height 6: 

,7.705 ,e2-0419 

^2 = e'^ ~ e"' 
Above, we wrote A^o as a tower of height 6: 

A^o ~ . 

The only difference is that in A'o, the top-level exponent is e instead of 2.0419. But e ~ 
2.0419 • 1.3313, so 

A^o ^ Sl-^^^\ 
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8. An Even Larger Number 



One can continue indefinitely playing this game: choose an arithmetic progression with an 
even more sparse distribution of primes; let be the sequence of reciprocals of those primes, 
then estimate the point at which the sum of these reciprocals exceeds 27r, so that equation (1.1) 
holds for more and more values of before it fails. 

With that in mind, here is our final example. 

As before, let P be the largest known prime, and let q = P^. We know that arithmetic 
progression qn + 1 contains infinitely many primes for n > 0. 

If p is any prime, then (f){p^) = — p^~^. Therefore, 0(g) = P^ — P^^^ = Q — P^^^, which 
proportionally speaking, is very close to q. Again, let's take M{q, 1) to be 0. 

The value of x for which the sum 



surpasses 27r is about 

X ~ exp(exp(0(g) ■ {2n — M{P, 1)) )) ~ exp(exp(27r(g — P^^^))) > exp(exp(7rg)). 

The inequality on the right holds for P > 2. Applying the same simplifying approximations in 
the previous section, the number of primes in this arithmetic progression up to this x is about 




p=i mod q 



X exp(exp(7rg)) exp(exp(7rg) — nq) 



0(g) log a; (g — P-^"^) exp(7rg) q — P^~^ 



Just as in the previous section. 




The fraction on the right is greater than 1. So, we will take as our approximation, 

A^o — exp(exp(g)) = exp(exp(P^)) 

where, approximately. 



pP ^ (loio^)io 



Just for fun, let's write this Nq as a tower of e's. We have 



logloglogA/'o ~ logP^ = PlogP ~ (10^°') ■ 10^ log 10 > (10^°') ■ 10^ 



It follows that 



log log log log log A^o > loglog(10^°' ■ 10^) ~ 16.95 > 



Therefore, we can write A^^o as, approximately. 



exp exp exp exp exp exp(e) = e' 



FUN WITH VERY LARGE NUMBERS 
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9. How Small is the Left Side Minus the Right Side of Equation (1.1)? 

When N is large enough that (1-1) does fail, the left side minus the right side is positive. 
As mentioned in Sections 2 and 4, Crandall has shown that this difference is surprisingly small 
when Ofc = l/{2k + 1), and when = prime. But how small are the differences for the 

ttfc that we consider in Sections 3, 5, 7, and 8? One would expect that the differences would be 
much smaller. The author has not pursued these questions yet. Estimating how much smaller 
they are could make an interesting exercise for the reader. 
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